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A study of the cold plasma cylindrical waves that may propagate in a specific type of 
two-dimensional magnetic field is initiated in this paper. The plasma is assumed to be of 
uniform density and collisionless, and a "spoke-wheel" magnetic field is considered which 
is both anisotropic and inhomogeneous (varying as the inverse radius), as defined in the text. 
Perturbation series solutions are obtained for the first Fourier component of the electric 
field for the four extreme cases: large and small magnetic field; large and small plasma 
densities. 

1. Introduction 

Interest in the analysis of wave propagation through magneto-active plasmas has, in the 
recent past, received renewed stimulation due to the emergence of the problem of communica- 
tion "black out" to or from reentry vehicles. 

With regard to this problem, in the present paper, we shall initiate a study of the cylindrical 
cold plasma [Allis, 1959; Auer et ah, 1958; Mason and Gold, 1962; Gold, 1963] waves that may 
propagate in a specific type of two-dimensional inhomogeneous magnetic field. The plasma 
is assumed to be of uniform density, and collisions are neglected. A so-called "spoke-wheel" 
magnetic field is considered which consists of two orthogonal (r, 6) components. The circular 

A m ... A 

component 2 is the field generated by a straight wire, while the remaining r component may 
may be approximated by the field generated at the end of a flattened solenoid of infinite radius. 

A A 

Both the r and fields decay as the inverse radius. The problems we shall investigate 

pertain to a specific (i.e., only r dependent, for the first Fourier component) class of waves 

which are sustained in a plasma that is situated in the region of two-dimensional symmetry of 

these fields. The related field lines are depicted in figure 1 (constant z). 

The r~ l dependence of the Be field generated by a straight wire is classical. The similar 

dependence of the B r field follows from simple geometrical consideration. If there are N B T 

lines issuing from the source of radius r , so that the B field there is aN/2Tr Q (where a is a 

constant), then at the radius r, the value of the field is aN!2wr, since the same number of lines 

cross all circles. 




Figure 1. Spoke-wheel magnetic field. 



i Assistant Professor of Physics, New York University, University Heights, New York, N.Y. 
2 The roof notation denotes a unit vector. 
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Similar studies of guided cylindrical waves in a cold plasma may be found in the literature 
[Astrom, 1950; Buchsbaum et al., I960]. However, for the most part, the steady B field in 
these investigations is in the ^-direction and is, of course, constant in space. The striking 
dissimilarity between such a B z field and the B T e fields herein stated is that in the former case 
the components of the related electric fields are decoupled, while in the latter case, even for 
constant-in-space B r6 fields, the E r j, z fields are severely coupled. A class of r-dependent 
solutions is examined for the magnetic field (I? r ,0,0). In this case, the E, field is completely 
decoupled from Ee, z and is nontrivial only if the signal frequency is equal to the plasma fre- 
quency. The (0,50,0) problem is not of current interest although it is amenable to a similar 
analysis. 

In this introductory discourse, we shall generate perturbation series for the electric fields. 
The two physical variables at our disposal are the plasma density and the steady magnetic 
field intensity. Of the four extremes that are naturally suggested, only the one of large plasma 
density prohibits perturbation analysis, since in that case the zero-order solution (infinite 
plasma density) yields zero electric field. In the remaining cases, quadratures are obtained 
for the first-order terms in the perturbation parameter. 

2. Analysis 

2.1. Equations 

The defining characteristic of a cold plasma is that the constitutive particles suffer no 
spread in velocity space, viz, the distribution function /(f) is a delta function <5({— £ n ). We 
shall consider also a two-component plasma whose ion mass far exceeds the electron mass so 
that in any subsequent perturbation from equilibrium, only the electronic motion contributes 
to the current J. The related wave equation for the electric field E (MKS units) is, 



(co 2 -o)^4- 1 -cVXrX)E=0 

where co is the mode frequency, co is the plasma frequency, and c is the speed of light. 
Let us consider the operator A in more detail. Two cases follow: 
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where we have set 
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If O is the magentic Larmor frequency [eB/m], and b is a unit vector B/B, then Kis the magentic 
wave number, defined through 
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In similar manner, there follows for A d 
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For the inverses we obtain 
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Returning to the general analysis, we may write (1) in component form. The depend- 
ence of the solution is absorbed into the circular harmonic form, 



E=S E'(r)e±'". 

1=0 

For case 1 (p=l/K T r), (1) appears as 

(^-o^-^E r -(±il) ^ A (rS#)=0, 



[(■^^vie^)]* 



"o£ 



l+^ : 



- a £,-c 2 (±M) 



•n <> //'- 



1(f)-. 



L r 1+F r or \ Or/ J 1+jr 



where E T j lias been written for E l r j. The equations constituting case 2 (j) = l/K e r) are similar 
in form but will not be considered in detail at t he present time. 

2.2. Solution 

In this analysis, w T e will restrict our considerations to the 1=0 mode. In this case, the 
E r field is completely decoupled from Ee, z and is nontrivial only for co=o? . The remaining 
two equations for E$ t z appear as : 



where we have set 



The L operators are 



so that 



[i+i z„(*)] e.-j^^e, 



a 2 =-^> y=Kr=p x ; x=k r; h=~* 



j 2 ^ I (f/ 2 



(2) 



(L n +x 2 )[AJ n (x)+BY n (x)}=0, 
where J n and Y u are nth order Bessel functions of the first and second kinds, respectively. 



327 



We will now consider four classes of perturbation solutions. The parameter of smallness 
for the four cases is defined in the following table : 

K 

(a) y=jr x =^ x ' Large B field (£ »i£) 

(b) y=§- x=e~ l x: Small B field (k «K) 

Kq 

(c) a 2 =e: Rare plasma (co «co) 

(d) a~ 2 =e: Dense plasma (co »co). 

Case a, large B field. In all of the following cases, the fields are expanded in powers of e : 

n=0 

In addition, all functions of e in the set (2) are expanded about e = 0. For the case under 
consideration, there results (y=ex), 

{x 2 +L 1 )^El n U n =a\ 2 x'(^]{-) n e 2n x 2n ) (j}El n >e n +j^ JZE[ n) e n \ 

(x 2 +L )^E^e n =a 2 e 2 x'(^(-) n e 2n x 2n ) (j}E™e n ~ J2E ( e n) e n \ 

The recursive] equations for [2? 2 (n) ; E { e n) ] are obtained by equating the coefficients of equal 
powers of e. The first two of these equations appear as (assuming that a is of order one) 

(x 2 +i 1 )^ (0 >=0; (x 2 +L )Ei*=0; (3a, b) 

{x 2 +L l )E^=a 2 x^E ( f ) ) {x 2 +U)E^ = -a 2 x z El°\ (4a, b) 

These two equations determine the solution to 0(e), which, in turn, appears as 

E=E (0) +eE (1 >; e=^=-^r; & = e ^; B ^r"\ 

The zeroth order solutions are the free space solutions, or, more accurately, the infinite 
magnetic field solutions. The simulation of the free space case is due to the rigid ordering 
of the particles by the B field, so that it becomes impossible for the electric fields to do work on 
them. For E (0) , we find, 

E2»=AJ (x)+BYo(x)=j ^ 

Er=CJ 1 (x)+DY 1 (x)=j 1 J 

where we recall x=k Q r, and J n and Y n are Bessel functions. 

The next order terms are composed of two parts : a particular integral and a complemen- 
tary function. The complementary functions for (4a) through (4b) are all the same as the so- 
lution to the zeroth order (3) . A particular solution of (4b) is obtained by setting 

There results the following first-order equation for<§ 2 ': 

&'+&^[ln(a^)]=-^ 
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where a prime denotes a derivative with respect to x. Multiplying through by exp [In (xjl)] = 
xjl gives 

5 (8 f z xfo) = -«Vml= - a ~ ± (xjj 2 , (6) 

the last equality following if the constants of jo are set equal to those of j x since (d/dx)(xji) =xj Q . 
Integrating (6) twice gives the desired result in quadrature form: 

In a similar manner, we obtain for Se (where E ( l ) = Seji), 

To order e the fields now appear as: 

E z =Ef(l + tg z ), 

E e =E^+eErge 

with the zero order solutions given by (5) and the functions §(x) given by (7). It should be 
noted first that the solutions to the first-order homogeneous equations have been absorbed in the 
unity factor, and second that the constants of E z m are not the same as those of Ef; 
however, A=C, B=D for the 0(e) solutions. Reference to (5) indicates that four arbitrary 
constants remain, two residing in E 2 and two in Eo. 

Case b, small B () field. This case is characterized by setting 

AuQ 

Equations (2) then appear as 



V + l o) e 1 =^[ tj ^](e,-*-e,) 



Expansion of the relevant functions about e=0 and ordering, results in a sequence of coupled 
equations for the field components E (n) - The first two equations appear as [for a = 0(1)1 

[L l + (l-a 2 )xW^ ) = 0;[L +(l-a 2 )x 2 }E^=0; 

[L 1 + (l-a 2 )x 2 ]El 1) =a 2 xE z 0) ;[L +(l-a 2 )x 2 }E 2 1) ^ 

The solutions to the zeroth order equations are given by 

Ei m =AJ 1 (a'x) + BY 1 (a'x) =j i; E™ = CMa'x)+DY (a'x) sj , 
where 

2 
/2 1 2 1 % 

a — 1 — a — 1 7/ 

or 

The next order solutions are obtained in a manner exactly the same as for case a; viz, we set 



There results: 



2a' '" 2a' Jx tfi(a't) 2a' \_ a 'tj^a 't)Jz 
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using the identity (d/dx)[j (a f x)]= —aj ] (a / x) and requiring that A= C, B=D. To within the 
stated order, the solution now appears as 

E 9 =E?(l+€g t ) 9 

E 9 =ES*+eE?g$. 

Case c, rare plasma. Here, we set 

In expanded form, (2) now appear as 

(z 2 +£0 2 E^e«=e (~^) (2 ^->«-+ir I 2 EJ->«-), 

For this case, it is particularly simple to write the general nth order recursive equations for E, 
which appear as [provided K/k = 0(1)]: 



(x 2 +L )El"™=Q^)(El»-y-Wp). 



In exactly the same manner as with the previous two cases, we obtain for E, to 0(e) 

E e =E^(l+ege), 

where 

Ep=j 1 (xy,E?>=j (x), 

and the first-order contributions are the quadratures, 

®*~J* *£(*) J*. 1+XV L ' X*Ji(*)J ' 

»._ f* <** p xvjgoo r. j,(g) "], 

62 J«,#)Ji, HXVL A«j>G0J ' 

where we have set X=K/k Q . 

Case d, dense plasma. This final case is characterized by the statment 

1 

ce- 
lt is easily seen that this equality renders the perturbation series 

E=E (0) +eE (1) + . . . 

degenerate if A is of order one, inasmuch as the starting term vanishes, generating, in turn, a 
null series. The physical significance is, of course, that electromagnetic waves cannot pass 
through an infinitely dense plasma (e=0). 
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3. Conclusions 

We have presented an introductory discussion on the class of electromagnetic waves 
which will pass through a specific type of cylindrical magnetic field embedded in a cold collision- 
less plasma of uniform density. The magnetic field is both anisotropic and inhomogeneous 
(varying as the inverse radius). The spoke-wheel magnetic field is defined in the text. 

The most significant aspect of the analysis is the presentation of a formalism (i.e., well- 
known Fourier expansion in the 6 variation together with subsequent perturbation expansion 
of the Fourier components) that permits the analytic investigation of these propagating fields. 
Four distinct cases were examined for the lowest order Fourier components in the four extremes: 
large and small magnetic field; large and small plasma densities. These in turn were examined 
to terms of first order in the parameter of smallness. The results are given in quadrature 
form. Well-known facts, such as infinite magnetic field being equivalent to vanishing plasma 
density, are recaptured in the lowest order expressions. The analysis differs from previous 
similar studies of guided electromagnetic weaves through cylindrical plasmas with embedded 
B z fields insofar as the electric field components in the present work are severely coupled while 
in the B z case they are not. 

It remains to fit, these solutions to specific geometries. The most physically relevant of 
these is the class of annular geometries with the steady B Q field zero in the region about the 
origin, and of spoke-wheel type exterior to this circular region. The plasma density is either 
zei'o or* nonzero (finite) in the various annular regions. 



The author wishes to express his deep appreciation to Bernard Wieder for his invaluable 
assistance rendered in the preparation of this paper. 

Much of the included analysis was completed under ( Contract AF-04(695)169 for Aerospace 
Corporation during the author's stay there in the summer of L962. 
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